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Abstr act. The extended Liiroth's Theorem says that if the transcendence de- 
gree of K(fi , . . . , f m )/K is 1 then there exists / 6 K(X) such that K(fi, . . . , f m ) 
is equal to K(/). In this paper we show how to compute / with a probabilistic 
algorithm. We also describe a probabilistic and a deterministic algorithm for 
the decomposition of multivariate rational functions. The probabilistic algo- 
rithms proposed in this paper are softly optimal when n is fixed and d tends 
to infinity. We also give an indecomposability test based on gcd computations 
and Newton's polytope. In the last section, we show that we get a polynomial 
time algorithm, with a minor modification in the exponential time decompo- 
sition algorithm proposed by Gutierez-Rubio-Sevilla in 2001. 



Introduction 

Polynomial decomposition is the problem of representing a given polynomial 
f(x) as a functional composition g(h(xj) of polynomials of smaller degree. This 
decomposition has been widely studied since 1922, see (22], and efficient algorithms 
are known in the univariate case, see [3] [201 133 EH] and in the multivariate case 

menus]. 

The decomposition of rational functions has also been studied, [HJ[I]. In the 
multivariate case the situation is the following: 

Let f(X x , ...,X n ) = h{X u . . .,X n )/f 2 {X 1} . . .,X n ) e K(X lt . .. ,X n ) be a ratio- 
nal function, where K is a field and n > 2. It is commonly said to be composite if 
it can be written / = u o h where h(Xi, . . . , X n ) £ K(Xi, . . . , X n ) and u G K(T) 
such that deg(u) > 2 (recall that the degree of a rational function is the maximum 
of the degrees of its numerator and denominator after reduction), otherwise / is 
said to be non-composite. 

This decomposition appears when we study the kernel of a derivation, see [2~i] . 
In [24] the author gives a multivariate rational function decomposition algorithm, 
but this algorithm is not optimal and works only for fields of characteristic zero. 
In this paper, we give a probabilistic optimal algorithm. In other words, our al- 
gorithm decomposes / € K(Xi, . . . ,X n ) with 0(d n ) arithmetic operations, where 
d is the degree of /. We suppose in this work that d tends to infinity and n is 
fixed. We use the classical O and O ("soft O") notation in the neighborhood of 
infinity as defined in [35] Chapter 25.7]. Informally speaking, "soft 0"s are used 
for readability in order to hide logarithmic factors in complexity estimates. Then, 
the size of the input and the number of arithmetic operations performed by our 
algorithm have the same order of magnitude. This is the reason why we call our 
algorithm "optimal" . 
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Furthermore, our algorithm also works if the characteristic of IK is greater than 
d(d- 1) + 1. 

This decomposition also appears when we study intermediate fields of an unira- 
tional field. In this situation, the problem is the following: we have m multivariate 
rational functions fi (X), . . . , f m (X) G ~K(X), and we want to know if there exists a 
proper intermediate field F such that K(fi, . . . , f m ) cFc KQQ. In the affirmative 
case, we want to compute F. If tr.deg K (F) = 1 then by the extended Liiroth's 
Theorem, see [23 Theorem 3 p. 15] we have F = K(/). 

Theorem 1 (Extended Liiroth's Theorem). Let ¥ be a field such that K C F C 
K(Xi, . . . , X n ) and tr.deg K (¥) = 1. Then there exists f £ K(Xi, . . . , X n ) such that 
F = K(/). 

The classical Liiroth's Theorem is stated with univariate rational functions. The- 
orem[T]gives an extension to multivariate rational functions. This extended theorem 
was first proved by Gordan in characteristic zero, see |13j . and by Igusa in general, 
see [T7J. There exist algorithms to compute /, called a Liiroth's generator, see e.g. 

USES. 

Thanks to the Extended Liiroth's Theorem the computation of intermediate fields 
is divided into two parts: first we compute a Liiroth's generator /, and second we 
decompose /. Then / = u o h, and F = K(h) is an intermediate field. In [15] the 
authors show that the decomposition of / bijectively corresponds to intermediate 
fields. They also give algorithms to compute a Liiroth's generator and to decompose 
it. Unfortunately, the decomposition algorithm has an exponential time complex- 
ity, but the complexity analysis of this algorithm is too pessimistic. Indeed, in the 
last section of this paper we show that we can modify it and get an algorithm with 
a polynomial time complexity. 

The decomposition of rational functions also appears when we study the spec- 
trum of a rational function. In this paper we use this point of view in order to give 
fast algorithms. 

Let K be an algebraic closure of K. Let / = /1//2 € K(Xi, . . . , X n ) be a rational 
function of degree d. The set 

h) = {(m : A) G Pjg; I M /i - A/a is reducible in K[Xl , . . . , X n ] , 

or degO/i - A/ 2 ) < d} 

is the spectrum of / = /i//2- We recall that a polynomial reducible in K[Xi , . . . , X n ) 
is said to be absolutely reducible. 

A classical theorem of Bertini and Krull, see Theorem [22l implies that u(/i, /a) is 
finite if /1//2 is non-composite. Actually, c(/i,/2) is finite if and only if /1//2 is 
non-composite and if and only if the pencil of algebraic curves /i/i — A/2 = has 
an irreducible general element (see for instance [HI Chapitre 2, Theoreme 3.4.6] 
and [7J Theorem 2.2] for detailed proofs). 

To the author's knowledge, the first effective result about the spectrum has been 
given bv Poincare [26], He showed that |cr(/i,/ 2 )| < {2d - l) 2 + 2d + 2. This bound 
was improved by Ruppert [28] who proved that 



\a(fx,h)\<d 2 -l. 
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This result was obtained as a byproduct of a very interesting technique developed 
to decide the reducibility of an algebraic plane curve. 
Several papers improve this result, see e.g. [2.3J 133 [21 El H] ■ 

The previous result says that if /1//2 is a non-composite reduced rational func- 
tion then for all but a finite number of A e K we have: f\ + A/2 is absolutely 
irreducible (i.e. irreducible in K[Xi, . . . ,X n }). Furthermore, the number of "bad" 
values of A is lower than d 2 — 1. Thus we can deduce a probabilistic test for 
the decomposition of a rational function, based on an absolute irreducibility test. 
In this paper we will give a decomposition algorithm based on this kind of idea. 
Furthermore, we will see that this algorithm is softly optimal when the following 
hypotheses are satisfied: 

Hypothesis (C): 
K is a perfect field of characteristic or at least d(d — 1) + 1. 

Hypothesis (H): 
deg(/i + A/ 2 ) = deg Xn (fi + A/ 2 ), where A is a new variable, 
{ii)Res Xn (fi(0,X n ) +A/ 2 (0,X n ), d x Ji(0,X n ) + A9 X „/ 2 (0,X„)) ^ in IK [A]. 

where deg x f represents the partial degree of / in the variable X n , deg / is the 
total degree of / and Resx n denotes the resultant relatively to the variable X n . 

These hypotheses are necessary, because we will use the factorization algorithms 
proposed in [22], where these kinds of hypotheses are needed. Actually, in [22] the 
author studies the factorization of a polynomial F and uses hypothesis (C) and 
hypothesis (L), where (L) is the following: 

Hypothesis (L): 

{(«) deg Xn F = deg F, and F is monic in X n , 
(ii)Res Xn (F(0,X n ),-^(0,X n )) ^ 0. 

If F is squarefree, then hypothesis (L) is not restrictive since it can be assured by 
means of a generic linear change of variables, but we will not discuss this question 
here (for a complete treatment in the bivariate case, see [TU] Proposition 1]). 

Roughly speaking, our hypothesis (H) is the hypothesis (L) applied to the poly- 
nomial /1 + A/2. In (H,i) we do not assume that f\ + A/2 is monic in X n . Indeed, 
after a generic linear change of coordinates, the leading coefficient relatively to X n 
can be written: a + Ab, with a, b £ K. In our probabilistic algorithm, we evaluate A 
to A ^ ct(/i, / 2 ), thus deg(/i + A/ 2 ) = deg(/i + A/ 2 ) and a + Xb ^ 0. Then we can 
consider the monic part of f\ + A/ 2 and we get a polynomial satisfying (L,i). Then 
(H,z) is sufficient in our situation. Furthermore, in this paper, we assume /i// 2 to 
be reduced, i.e. f\ and fi are coprime. We recall in Lemma[5]that in this situation 
fx + A/2 is squarefree. Thus hypothesis (H) is not restrictive. 

Complexity model. In this paper the complexity estimates charge a constant 
cost for each arithmetic operation (+, — , x, -=-) and the equality test. All the 
constants in the base fields (or rings) are thought to be freely at our disposal. 

In this paper we suppose that the number of variables n is fixed and that the 
degree d tends to infinity. Furthermore, we say that an algorithm is softly optimal 
if it works with O(N) arithmetic operations where N is the size of the input. 
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Polynomials are represented by dense vectors of their coefficients in the usual 
monomial basis. For each integer d, we assume that we are given a computation tree 
that computes the product of two univariate polynomials of degree at most d with 
at most 0(d) operations, independently of the base ring, see [39] Theorem 8.23]. 
We use the constant uj to denote a feasible matrix multiplication exponent as defined 
in [3§] Chapter 12]: two n x n matrices over K can be multiplied with C(n") field 
operations. As in [5] we require that 2 < uj < 2.376. We recall that the computation 
of a solution basis of a linear system with m equations and d < m unknowns over 
K takes ©(mcT" 1 ) operations in K Chapter 2] (see also [33 Theorem 2.10]). 
In [35] the author gives a probabilistic (resp. deterministic) algorithm for the 
multivariate rational factorization. The rational factorization of a polynomial / 
is the factorization in K[X], where K is the coefficient field of /. This algorithm 
uses one factorization of a univariate polynomial of degree d and 0{d n ) (resp. 
(5(ci ra+w-1 )) arithmetic operations, where d is the total degree of the polynomial 
and n > 3 is the number of variables. If n = 2, in [21], [HI Errata], the author gives 
a probabilistic (resp. deterministic) algorithm for the rational factorization. The 
number of arithmetic operations of this algorithm belongs to <D(d 3 ) (resp. 0(d u+1 )). 
We note that for n > 3 if the cost of the univariate polynomial factorization belongs 
to 0(d n ) then the probabilistic algorithm is softly optimal. 

Main Theorems. The following theorems give the complexity results about our 
algorithms. Although we will use no probabilistic model of computation, we will 
informally say probabilistic algorithms when speaking about the computation trees 
occurring in the next theorems. For the sake of precision, we prefer to express the 
probabilistic aspects in terms of families of computation trees. Almost all the trees 
of a family are expected to be executable on a given input (if the cardinality of IK 
is large enough). 

Theorem 2. Let f = fi/fs be a multivariate rational function in K(A"i, . . . , X n ) 
of degree d, there exists a family of computation trees over K parametrized by z :— 
(a, b) G K 2n such that: 

• Any executable tree of the family returns a decomposition u o h of f with h 
a non-composite rational function. 

• If®ik are n °t the roots of some non-zero polynomials the tree corresponding 
to z is executable. 

Furthermore, we have: 

(1) An executable tree performs two factorizations in K[Ai, . . . ,X n ] of polyno- 
mials with degree d, and one computation of u. 

(2) Under hypothesis (C) and (H) we have this estimate: an executable tree 
performs one factorization of a univariate polynomial of degree d over K 
plus a number of operations in K belonging to 0(d n ) if n > 3, or to 0(d 3 ) 
ifn = 2. 

Since we use the dense representation of f\ and fi , the size of / is of the order 
of magnitude d n . The previous statement thus asserts that the complexity of our 
probabilistic algorithm is softly optimal for n > 3. 

We precise the condition "If a, b are not the roots of some non-zero polynomials" 
in Remark [13] and Remark [15] 

In characteristic zero we can say that for almost all z the tree corresponding to z 
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is executable. 

We also give a deterministic decomposition algorithm. 

Theorem 3. IfK is a field with a least max(d 2 , |d 2 — 2d + 1) elements, then the 
decomposition f = uoh, with h non- composite, can be computed with at most 0(d 2 ) 
absolute factorizations of polynomials with degree d, and at most 0(d 2 ) computa- 
tions of u where f and h are given. 

If we can use the algorithm proposed in [TU] and [22 , as we will see in Remark 
II 81 our deterministic algorithm uses one factorization of a univariate polynomial of 
degree d with algebraic coefficients of degree at most d, and at most 0{d n+u+2 ) if 
n > 3 or <D{d 6 ) if n = 2 arithmetic operations in K. 

With the tools used for the decomposition algorithms, we can compute a Liiroth's 
generator. 

Theorem 4. Let fi, . . . ,f m G K(Ai, . . . , X n ) be m rational functions of degree 
at most d. There exists a family of computation trees over IK parametrized by 
z = (zi, . . . , z m ) G K 2nm , such that: 

If for all i = 1, . . . ,m, zi E K 2 ™ belongs to an open Zariski set related to ft, . . . ,fj 
then the tree corresponding to z is executable on ft, ... ,f m and it returns a Luroth 's 
generator of K(f i , . . . , f m ) . 
Furthermore, we have: 

(1) An executable tree performs 2m gcd computations in K[-Xi, . . . ,X n ] with 
polynomials of degree at most d. 

(2) 7/K has at least (Ad+2)d elements then we have the estimate: an executable 
tree performs 0(md n ) arithmetic operations in IK. 

As before, this algorithm is softly optimal because the order of magnitude of the 
input is tog?™. A precise description of the open Zariski set is given in Remark 1291 

In the last section we prove the following result: 

Theorem 5. Let f = f x /f 2 G K(X). 

f = uoh, with h = ht/ha if and only if H(X_,Y) = ht(X_)h 2 (Y) - h 2 (2Qh 1 (Y) 
divides F{X,Y) = h(X)f 2 (Y) - f 2 (X)f 1 (Y). 

Furthermore, if h\jh 2 is a reduced non-composite rational function then hi is one 
of the irreducible factors with the smallest degree relatively to X_ of F. 

The first part of this theorem is already known, see [3D] ■ Here, we prove that H 
is irreducible if h\jh 2 is non-composite. This result implies that we can modify the 
exponential time decomposition algorithm presented in |15) and get a polynomial 
time algorithm. 

Comparison with other algorithms. There already exist several algorithms for 
the decomposition of rational functions. In |15j . the authors provide two algorithms 
to decompose a multivariate rational function. These algorithms run in exponential 
time in the worst case . In the first one we have to factorize ft {20 f 2 QQ — ft (JQ / 2 {20 
and to look for factors of the following kind hi (X)h?. (Y) — hi {Y)h 2 {20 ■ The authors 
say that in the worst case the number of candidates to be tested is exponential in 
d = deg{fi/f 2 ). In the last section we show that actually the number of candidates 
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is bounded by d. Thus we can get a polynomial time algorithm. 
In the second algorithm, for each pair of factors (/i.1,/12) of /1 and fi (i.e. hi di- 
vides /1 and h 2 divides /2), we have to test if there exists u £ K(T) such that 
/1//2 = u(h\/h,2)- Thus in the worst case we also have an exponential number of 
candidates to be tested. 

To the author's knowledge, the first polynomial time algorithm is due to J.-M. 01- 
lagnier, see [24] ■ This algorithm relies on the study of the kernel of the following 
derivation: S U (F) = uo A dF, where F £ K\X] and w = f2d.f1 - fidf 2 . In [53] 
the author shows that we can reduce the decomposition of a rational function to 
linear algebra. The bottleneck of this algorithm is the computation of the kernel 
of a matrix. The size of this matrix is 0(d n ) x 0(d n ), then the complexity of this 
deterministic algorithm belongs to 0(d nuj ). In [24], as in this paper, the study of 
the pencil — A/2 plays a crucial role. 

Structure of this paper. In Section 1, we give a toolbox where we recall some 
results about decomposition and factorization. In Section 2, we describe our al- 
gorithms to decompose multivariate rational functions. In Section 3, we give an 
indecomposability test based on the study of a Newton's polytope. In Section 4, 
we give two algorithms to compute a Liiroth's generator. In Section 5 we show that 
the decomposition algorithm presented in |15j can be modified to get a polynomial 
time complexity algorithm. 

Notations. All the rational functions are supposed to be reduced. 
Given a polynomial /, deg(/) denotes its total degree. 
K is an algebraic closure of K. 

For the sake of simplicity, sometimes we write K.[X] instead of K[A"i, . . . ,X n ], for 
n > 2. 

Res(A, B) denotes the resultant of two univariate polynomials A and B. 
For any polynomial P £ K[X], we write U{P) := {a £ K n | P(a) ^ 0}. 

1. Prerequisite 

The following result implies, as mentioned in the introduction, that hypothesis 
(H) is not restrictive. 

Lemma 6. J//1//2 is reduced in K(Xi, . . . ,X n ), where n>l and A is a variable, 
then fi + A/2 is squarefree. 

Now we introduce our main tools. 

Proposition 7. Let f = /1//2 be a rational function in K(Ai, . . . , X n ). 

f is composite if and only if nf\ — A/2 is reducible in K[X] for all fj,, X £ K such 

that deg(/x/i - A/ 2 ) = deg(/). 

We also have: f is non-composite if and only if its spectrum cr(/i, $2) is finite, 
if and only if fi — T/2 is absolutely irreducible in K.(T)[X], where T is a new 
variable. 

Furthermore if deg(/) = d then c(fi,f 2 ) contains at most d 2 — 1 elements. 

Proof. The first part of this result was known by Poincare see [26] , for a modern 
statement and a proof, see [7] Corollary 2.3]. 

The bound | a (hi, | < d 2 — 1 is proved for any field in the bivariate case in 
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[23j . We deduce the multivariate case easily thanks to the Bertini's irreducibility 
theorem, see e.g. [7] or the proof of Theorem 13 in [4] for an application of the 
Bertini's irreducibility theorem in this context. 

□ 

Lemma 8. Let h = hi/h 2 be a rational function in KQQ, u = ui/u 2 a rational 
function in K(T) and set f = u a h with f = /1//2 G 'K(X). For all A G K such 
that deg(ui — XU2) — degu, we have 

fi - A/ 2 = e(hi - tih 2 ) ■■■(hi- t k h 2 ) 

where e £ K, k — degu and ^ G K are the roots of the univariate polynomial 
ui(T)-Au 2 (T). 

Proof. See the proof of Lemma [39] in Section [5] Lemma [39] is a generalization of 
Lemma [8] We state Lemma [8] in our toolbox because the generalization will be 
only used in Section □ 

Remark 9. If ti € K then hi — tih 2 <E K[-Xi, . . . , X n ] is an irreducible factor of f\ — 
A/2. Thus with a rational factorization we get information about the decomposition 
of /. This remark will be used during our probabilistic decomposition algorithm in 
order to avoid an absolute factorization. 



2. Decomposition algorithms 

2.1. Computation of u. Suppose that / = /1//2 = u o h e K(Xl, . . . , X n ), 
h e K(Xi, . . .,X n ), and u 6 K(T). We set h = hi/h 2 . 

Usually, when hi and h 2 are given we get u — ui/u 2 by solving a linear system, see 
[T51 Corollary 2]. Let M(h\, h 2 ) be the matrix corresponding to this linear system 
in the monomial basis. In our situation the size of M.(h\, h 2 ) is 0(d n ) x 0(d). Thus 
we can find u with C(d n+ " _1 ) operations in K. 

We can get u with another approach. This approach is based on a strategy due to 
Zippel in [IT] . Zippel showed in the univariate case that we can compute u quickly. 
His strategy is the following: compute the power series H such that hoH(X) = X, 
then compute foH, and finally deduce u with a Pade approximant. All these steps 
can be done with 0(d) or 0(d 3 / 2 ) arithmetic operations, see [HI Chapter 1], and 
[6]. Thus we deduce that in the univariate case, u can be computed with 0(d 3 ^ 2 ) 
arithmetic operations. 

In the multivariate case with hypothesis (H), we have deg(/) = dcg Xn (f). Thus 
f(0,X n ) = uo h(0, X n ) is not a constant. Then we can apply Zippel's strategy to 
/(0, X n ) in order to find u. This method is correct because if / and h are given 
then there is a unique u such that / = u o h, see [151 Corollary 2]. Thus we have 
proved the following result: 

Lemma 10. Let f,h £ K.(Xi, . . . , X n ) be rational functions. We suppose that f 
satisfies hypothesis (H) and we set deg( f) = d. If there exists u € K(T) such that 
f = u o h then we can compute u with 0(d n ) arithmetic operations. 

Proof. We compute /(0, X n ) with 0(d n ) arithmetic operations. Then we compute 
u as explained above with 0(d 3 ^ 2 ) arithmetic operations. □ 
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2.2. A probabilistic algorithm. Decomp 

Input: / = hjh G K(Xi, . . . ,X„), z := (a,b) e K 2n . 

Output: A decomposition of / if it exists, with / = u o h, u — Ui/112, h — hi/h-z 
non-composite and degu > 2. 

(1) We set F a = / 2 (a)/i(X) - h{a)f 2 {X), F b = f 2 {b)h{X) - fi(b)fc(X). 

(2) Factorize F a and Fb. 

(3) If F a or Fb is irreducible then Return "r is non-composite" . 

(4) Let T a (resp. J-b) be an irreducible factor of F a (resp. Fb) with the smallest 
degree. 

(5) Set h = TalTb- 

(6) Compute u such that / = u o h as explained in Section [2. II 

(7) Return u, /i. 

Exemple 11. 

a- We consider / = / x // 2 , with /1 = X 3 + Y 3 + 1 and f 2 = 3XY. We set 
a= (0,0), 6= (0,1). Then F a = and F 6 = 3X 3 + 3r 3 -6XF+3. F Q 

is reducible but F b is irreducible then we conclude that / is non-composite. 

b- Now, we apply the algorithm Decomp to the rational function / = u o h, 
where u = (T 2 + l)/T and h = h 1 /h 2 with hi = X 3 + Y 3 + l and h 2 = 3XY. 
We have seen above that h is non-composite. 
In this situation with a = (0, 0) and b = (0, 1) we get: 

F a = -3.X.Y.(X 3 + Y 3 + 1), and 

F b = -12.X.Y.(X 3 + Y 3 + 1). 

Then the algorithm cannot give a correct output in this situation. Here, 
we have /2(a) = f2{b), we will see that we must avoid this situation. 
If we set a = (2, 1) and b = (1, —1) then: 

F a = 60. {X 3 + Y 3 - 5XY + 1).(X 3 + Y 3 - -XY + 1), and 

5 

F b = -3.{X 3 + Y 3 + XY + l).(X 3 + Y 3 + 3XY + 1). 

Then we get J- a — X 3 + Y 3 — 5XY + 1 and T h = X 3 + Y 3 + XY + 1. 
The algorithm Decomp returns h = T a jTb- This is a correct output since 
UoF a /F b = hi/hn, where U = (T/6 + 5/6)/( - T/2 + 1/2). 

Proposition 12. If a,b are not the roots of some non-zero polynomials then the 
algorithm corresponding to z — (a, b) is correct. 

Proof. First, we suppose that / is non-composite and we set 

Spect A>/2 (T 1 ,T 2 )= J] (pT 2 - XTx). 

( M :A)G<t(/ 1 ,/ 2 ) 

We have Spectj * (^, A) = if and only if (/i : A) G o-(fi,f 2 ). 

IfSpect /i /2 (/ 2 (a),/i(a)).Spect /l /2 (/ 2 ©,/i(6)) ^ 0thenF a and F b are absolutely 

irreducible and degF a = degFb = deg/. 

This gives: if a and b avoid the roots of 

S(A,B) := S P ect /ii/2 (/ 2 (A),/ 1 U)).Spect /ii/2 (/ 2 (F),/ 1 (S)), 
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where degS* < 2d(cP — 1) by Proposition [3 then the algorithm returns: "r is non- 
composite" . 

Second, we suppose / = voH , with H G K(Xi, . . . , X n ) a non-composite rational 
function. We set v = V\ jv 2 , H = H\/ H 2 such that these two rational functions are 
reduced. We also suppose that /2(a) and f 2 (b) are nonzero. 

If deg F a — deg Ft, — deg / then a and b are not the roots of a polynomial D of 
degree d. Thanks to Lemma [S] we have: 

F a = e{H x - t x H 2 ) ■■■(H 1 - t k H 2 ), 

F b = e'{H x - s x H 2 ) -..(Hi- s k H 2 ), 

with e, e' S K, ij, Sj € K. 

As H 1 {a)/H 2 (a) (resp. H 1 {b)/H 2 {b)) is a root of f 2 (a)v 1 {T) - h{a)v 2 {T) (resp. 
Hk)vi{T) - fi(b)v 2 {T)), we set ti = Hi(a)/H 2 (a) and si = Hi(b)/H 2 {b), and we 
remark that ti, si £ K. We set 

S P ect ffliff2 (T)= J] (T-A). 

Aecr(Hi,H 2 )nK 

If Spect Hl t H 2 (h) ^ (resp. Spect ffl , H2 (si) ^ 0) then i?i - t\H 2 (resp. i/i —s\H 2 ) 

is absolutely irreducible. 

If 

i?(a,6) = Res T (f 2 {a) Vl {T) - h(a)v 2 (T), f 2 (b)v x {T) - fi(b)v 2 (T)) ? 

then ti ^ Sj for all i, j. We remark that i? is a nonzero polynomial by Lemma [6] 
since V\ and v 2 are coprime. Thus step H] gives J- a = Hi — tH 2 , T b = H\ — sH 2 
with t, s G K and t =/= s. Then /i = Ta/J~b is non-composite, because H\/H 2 is 
non-composite. □ 

Remark 13. Now, with the notations of the previous proof, we can explain in 
details the meaning of: "If a,b are not the roots of some non-zero polynomials" in 
Proposition [T^l and Theorem [21 This means: 
If / is non-composite then there exists a nonzero polynomial 

P(£,R) := S(A,B) 

of degree at most 2d(d 2 — 1) such that for any (a, 6) £ U{P) the algorithm corre- 
sponding to z is executable and returns a correct output. 
If / is composite then there exists a nonzero polynomial 

Dx{A,B) := f 2 (A). f 2 (B) . D (A) . D (B) 

of degree at most Ad such that; 

for any (a, b) € IA(D\), there exist nonzero polynomials 

D 2 (A);= 11 (H 2 (A) - XHx(A)) 

\ea(H 1 ,H 2 )r\K 

of degree at most (d 2 — l).d/2, and 

R(A,B) 

where deg^i? < d 2 /2 and deg_g R < d 2 /2, such that; for any (a, b) e U(D 2 {A).D 2 (B).R(A, B)) , 
the algorithm corresponding to z = (a, b) is executable and returns a correct output. 
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Proposition 14. Under hypotheses (C) and (H), if a and b are not the roots of 
a non-zero polynomial then we can use the algorithm proposed in [22j . Then the 
algorithm Decomp performs one factorization of a univariate polynomial of degree 
d over IK plus a number of operations in K belonging to 0{d n ) if n > 3 or to 0(d 3 ) 
ifn = 2. 

Proof. As / satisfies (H,z), we deduce that if a and b are not the roots of a polyno- 
mial D of degree d, then the monic part relatively to X n of F a (resp. Fb) satisfies 
(L,t). 
We set: 

P(A) = Resx n (h{0,X n ) - A/ 2 (0, A„), d Xn fi{0,X n ) - Afe n / 2 (0,X„)). 

By hypothesis (H,m), P(A) ^ in IK [A]. Furthermore if /2(a) and f 2 (b) are nonzero 
and ^(/i (a) //2(a)) 7^ (resp. f(/i®//2®) 7^ 0) then hypothesis (L,m) is satis- 
fied for F a (resp. Then we can use Lecerf's algorithm, see 22 . This gives: if 
a and b avoid the roots of 

V(A,B) = V(f 1 (A)/f 2 (A)).V(f 1 (B)/f 2 (B)).(f 2 (A)J 2 (B)) dcsV+ \ 

and degV < 2{d{d-l)d + d) then we can use the algorithm proposed by G. Lecerf 
in [22]. 

The complexity result comes from Lemma [TQl and f22j Proposition 5], [21] Propo- 
sition 2] and [22j Errata]. □ 

Remark 15. The meaning of the condition "if a and b are not the roots of a non- 
zero polynomial" in Proposition Q3] is the following: If we want to use Lecerf's 
factorization algorithm in order to get the complexity estimate given in the second 
part of Theorem [21 then a and b must also avoid the roots of the polynomial 

D(A).D(B).V(A,B), 
where deg£> < d and dcgV < 2(d 2 (d - 1) + d). 

It follows that Theorem [2] comes from Proposition [T2l and Proposition [T4l 
2.3. A deterministic algorithm. Decomp Det 

Input: / = /1//2 £ K.(Xi, . . . , X n ), S = {s , . . . , sg} a subset of K with at least 
B + 1 = max(d 2 , |d 2 — 2d + 1) distinct elements. 

Output: A decomposition of / if it exists, with / = u o h, u — ui/u 2 , h = h\jh 2 
non-composite and degu > 2. 

t:=false, A := 0. 
While t=false do 

(1) If deg(/i + s x f 2 ) = deg(/) then go to step [2] else A := A + 1. 

(2) Compute the absolute factorization of Fx ■= f\ + s\f 2 . 

(3) If F\ is absolutely irreducible then Return "/ is non-composite" . 

(4) If F\ is absolutely reducible then 

(a) If two distinct absolute irreducible factors fi, f 2 belong to K\X] then 
we set h\ := fi and h 2 := f 2 , 

If there exists an absolute irreducible factor fi := T\ + eF 2 , with 
££l\K and F\,J- 2 G K[X] then we set hi := Ti, h 2 := F 2 , 
Else A := A + 1 and go to stepQ] 



DECOMPOSITION OF MULTIVARIATE RATIONAL FUNCTIONS 



11 



(b) Compute u (if it exists) such that / = u o h as explained in Section 

ED 

(c) If u exists then t:=true else A := A + 1. 
Return u,h. 

Exemple 16. 

a- We consider / = f x /f 2 , where fx = 3XY and f 2 = X s + Y 3 + 1. This 
gives Fq = 3.X.Y, then Fo is reducible, and this gives h = X/Y . We do 
not find a rational function u such that / = u o (X/Y) then we consider 
F% = fi + f 2 . F% is absolutely irreducible, then the algorithm Decomp Det 
returns / is non-composite. 

b- Now, we apply the algorithm Decomp Det to the rational function / = uoh, 
where u — (T 2 + 1)/T and h — (X 3 + Y 3 + l)/(3XY). As we have seen 
above h is non-composite. 
In this situation we have: 

F = (X 3 + Y 3 + 1 + 3.i.X.Y)(X 3 + Y 3 + 1 - 3.1.1.7), 

where i 2 = — 1. 

Then we have fi = X 3 + Y 3 + 1 + 3.LX.Y, Fx = X 3 + Y 3 + 1, I 2 = 3XY. 
The algorithm returns J-1/J-2 = h. 

Proposition 17. The algorithm is correct. Furthermore we go back to step\l\at 
most 0(d ) times. 

Proof. First, we suppose that / is non-composite. By Proposition [7] there exists 
s\ £ S such that s\ £ o-(fx,fc) because S contains at least d 2 elements. Thus 
fx + S \ Q H is absolutely irreducible and step [3] returns / non-composite. 
We remark that if fx + s\f 2 is reducible then we cannot find u during step l4bl 
because / is non-composite. Then if / is non-composite the algorithm is correct. 

Second, we suppose that / is composite and / — v o H with H = Hx/Hz a 
reduced and non-composite rational function, deg v > 2 and v = Vx jvi is a reduced 
rational function. 

fx + s\f 2 = el\i{Hx +UH 2 ) by LemrnaH where (vx + sxV 2 ){ti) = 0. 
There exists s\ € S such that D(s\ ) 7^ 0, where 

D(A) — Res(vx + Av 2 , v[ + Av' 2 ) x Jj (v 2 (xi) - Avxfa)) . 

x i e<T(H 1 ,H 2 )nK 

Indeed D( A) is a nonzero polynomial by Lemma |6] since vx and v 2 are coprime. 
Furthermore, by Proposition we have 

degL* < degu(degw — 1) + ((degi/) 2 — l). deg v. 

As deg v. deg H = d and deg v > 2 , we get 

degD < 3/2d 2 - 2d. 

As S contains at least 3/ 2d 2 — 2d + 1 distinct elements, there exists s\ £ S such 
that D(s\ a ) 7^ and then for all i, t\ cr(Hx, H 2 ), and U 7^ tj for all i 7^ j. 
Then for Ao we construct hx and h 2 as explained in step 0a] (If tx,t 2 S K are 
distinct then we have two absolutely irreducible factors in IK[A], else if tx 6l\I 
then we construct hx and h 2 with only one absolutely irreducible factor.) We have 
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hi/fi2 = w o H1/H2 where w E K(T) and deg w = 1. 

We remark that if / is composite then we find a decomposition f = u o h with h 
non-composite. Indeed, there exist (fi : A) and (// : A') ^ (// : A) E Pj^ such that 
/i/ii + A/12 and /i'/ii + A'/i2 are absolutely irreducible. (It is obvious if tx, t% E K. If 
ii E K\K there exists a conjugate t[ of ^! over K such that hi + t[h 2 is absolutely 
irreducible.) Then hi/h 2 is non composite by Proposition [7] Thus if / is non- 
composite the output is correct. □ 

Theorem [3] is a direct corollary of Proposition [T7l 

Remark 18. In |10j the authors show that we can compute, under the hypothesis 
(C), the absolute factorization of a bivariate squarefree polynomial with at most 
0(d A ) arithmetic operations. As we go back to step 1 at most O(cP) times we 
deduce that the algorithm Decomp Det uses at most 0(d 6 ) arithmetic operations. 
When n > 3, a complexity analysis of an absolute factorization algorithm as studied 
in [10] is not done, but we can estimate the cost of our deterministic algorithm. In- 
deed, we can reduce absolute factorization to factorization over a suitable algebraic 
extension K[a] of degree at most d over K, [3H [3H [TH]. With this strategy and 
with the deterministic factorization algorithm proposed in [22] we get an absolute 
factorization algorithm which performs at most CD(d n+ul ~ 1 ) arithmetic operations 
in K[a}. Thus the algorithm performs 0(d n+u) ) arithmetic operations in K, because 
[K[a] :K]<d. As we go back to step 1 at most 0(d 2 ) times we deduce that, if we 
can use Lecerf's deterministic factorization algorithm, the algorithm Decomp Det 
uses at most 0{d n+ul+2 ) arithmetic operations and one factorization of a univariate 
polynomial of degree d with coefficients in K[a]. 



3. An indecomposability test using Newton's polytope 

In Section [21 if fx and f% are sparse our algorithms do not use this information. 
In this section we give an indecomposability test based on some properties of the 
Newton's polytope. The idea is to generalize this remark: if deg / is a prime 
integer then / is non-composite. This is obvious because f = uo h implies deg / = 
deg u. deg h, and deg u > 2. 

Definition 19. Let f(X) € K[X U . . .,X n ], the support of f(X) is the set S f of 
integer points . . . , i n ) such that the monomial AJ 1 ■ • • A^" appears in / with a 
nonzero coefficient. 

We denote by N(f) the convex hull (in the real space R") of Sf. This set N(f) is 
called the Newton's polytope of /. 

Definition 20. We set N(f\/ ' fi) = N(fi — A/2) where A is a variable, and where 
fx — A/2 is considered as a polynomial with coefficients in IK [A]. 

Remark 21. As A is a variable N(fx — A/2) is the convex hull of S/ 1 U Sf 2 . 

We recall the classical Bertini-Krull's theorem in our context, see [311 Theorem 
37]. 

Theorem 22. (Bertini-Krull) Let fx/ fi a reduced rational function. Then the 
following conditions are equivalent: 
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(1) /1//2 is composite, 

(2) (a) either there exist hi , h 2 G K[2d with deg y A (X)— A fy.(X) > max(deg /ii, deg ft. 2 ) 

and Oi(A) G K[A], smc/i </ia< 

e 

A (20 - a/ 2 Q0 = X^^^QO 6 ^; 

i=0 

(b) or the characteristic pofWLis positive and A (X)—A f?.(X) G K[A] [Xf , . . . , X£] . 

Lemma 23. If fx/ f 2 is a composite rational function and the characteristic p of 
K zs such that p = or p > d, then there exist e G N, /ii,/i 2 € K[X] swc/i i/iai 
N(fx/fc) = eN(hx/h 2 ). 

Proof. By Theorem H we have /i(X) -A/ 2 (X) = £* =0 a;(A)/ii(X) l fr 2 (X) e -\ We 
denote by it(A, x) the polynomial 

e e 

tt(A, X ) = e a »( A )^ = a ^ A ) n u - ^( a )) > 

i=0 i=l 

where <£j(A) G K(A). 
Thus 

e 

ji(X) - A/ 2 (X) - 0e (A) J] (/11(A) - <£,(A)/i 2 (X)) • 

i=l 

All the factors hi(X_) — (A:(A)/t?(X) G K(A)[X] have the same support. 

Indeed, if we suppose the converse then there exist a coefficient c\ G K of h\ and a 

coefficient c 2 G K of h 2 and two indices i and j such that: 

ci — ¥>i(A)c 2 = 0, ci — <£j(A)c 2 ^ 0. 

Then c 2 ^ and y>*(A) = Ci/c 2 G K. Thus /ii - tpi(A)h 2 G K[X] is a factor of 
fi{X) - A/ 2 (X). This implies /i(X) - A/ 2 (X) is reducible in K[A][X]. This is 
impossible because f\ and / 2 are coprime. 
Then, for alH = 1, . . . , e, we have: 

N(hi - ^(A)h 2 ) = N{hi - Ah 2 ) = Nfa/hi). 

We recall that F = F 1 .F 2 implies N(F) = N(Fx) + iV(F 2 ), see for example O 
Lemma 5], where the sum is the Minkowski's sum of convex sets. Thus we have: 

e e 

N{fi/f 2 ) = N(h - Af 2 ) = Y J N{h 1 - <pi(A)h 2 ) =Y,N(hM = eN{h x /h 2 ). 

i=l »=1 

This is the desired result. □ 

The previous lemma says that if / is composite then all the vertices of N(f) 
have a common factor: e. This gives our indecomposability test designed for sparse 
polynomials fx and f 2 : 

Corollary 24 (Indecomposability test). Letp be the characteristic ofK., andp = 
or p > d. 

Let (ix , • ■ • 7 4i )> • ■ • , (i[ k \ • • ■ 7 i^) be the vertices of N(f). 
//gcd(i^ 1 \ . . . , in \ ■ • ■ , i[ k \ • ■ • , in ) = 1 then f is non-composite. 
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4. Computation of a Luroth's generator 

In this section we show how to compute a Luroth's generator. We give two 
algorithms. The first one follows the strategy proposed in [35] for univariate rational 
functions. The second one uses the algorithm Decomp and the computation of a 
greatest common right component of a univariate rational function. 

4.1. Generalization of Sederberg's algorithm. In this subsection, we gener- 
alize Sederberg's algorithm. Sederberg's algorithm, see [32], is a probabilistic al- 
gorithm to compute a Luroth's generator in the univariate case. Here, we show 
that the same strategy works in the multivariate case. Our algorithm is also a 
kind of probabilistic version of the algorithm presented in [15] . Indeed, here we 
compute gcd of polynomials of the following kind /2(fi)/i(2Q — /i (0/2(20, where 
a £ K n . In [15], the authors compute gcd of polynomials of the following kind 
/2Q3/1QQ — /1 (H)./2(20, where Y_ are new independent variables. 

Sederberg Generalized 

Input: f(X) = /i// 2 (20, 3(20 = SiAte(20 S K(X U ...,X n ) two reduced rational 
functions, a, b £ K n , n > 2 . 

Output: h(X) £ K(X_) such that K(f,g) = K(h), if h exists. 

(1) F a := fo{a)h{X) - A(a)/ 2 (20, G a := g 2 (a) 9l (X) - 9l (a)g 2 (X). 
H a :=gcd(F a ,G a ). 

If H a is constant then Return "No Luroth's generator", else go to [5] 

(2) F b := h{b)h{X) - /i(6)/ 2 (20, G b := g 2 (b)g 1 (X) - 9l (b)g 2 (X). 
H b :=gcd(F b ,G b ). 

If H b is constant then Return "No Luroth's generator" , else go to [3J 

(3) Return h := H a /H b . 

Exemple 25. a- We set / = X, and g = Y, a = (0,0), b = (1,0). Thus 
F a = X, G a = Y and H a — 1. The algorithm Sederberg Generalized gives 
K(f, g) = K(X, Y) has "No Luroth's generator" . 
b- We consider f = Uoh and g = V oh where h = (X 3 + Y 3 + 1)/(3XY), 
U = T 2 /(T + 1), V = (T + 2)/(T 3 + 3). /lis a non-composite rational 
function. 

We set a = (0, 0), b = (2, 1). In this situation we have: 

H a = 3XY, and H b = 12. (X 3 + Y 3 - hXY + 1). 

The algorithm Sederberg Generalized returns H a /H b . This is a correct out- 
put because K(f,g) = K(h) and h = u o (H a /H b ) where u is the rational 
function u = (20T + 1)/(12T). 

Now, if we set a = (0, 0), b = (0, 1) then we get H a = 3XY and H b = 12XY. 
In this situation the output H a /H b is not correct. We are in a situation 
where h{a) = h(b) and we will see that we must avoid this situation. 

Proposition 26. There exists an open Zariski set U C K 2n related to fi and f 2 , 
such that for all (a, b) £ U the tree corresponding to (a, b) is executable on f, g and 
returns (if it exists) h such that K(/i) = K(/, g). 

In order to prove this proposition we recall some results. 
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Definition 27. Given fi,.. . ,f m G K(20, we sav that they have a common right 
component (CRC) h, if there are rational functions Ui G K(T), i = 1, . . . , m, such 
that fj = Ui o /i, and degiti > 1. 

/i is a greatest common right component (GCRC) of fi, . . . , f m if the m^s have not 
a common right component of degree greater than one. 

Proposition 28. K(fi, . . . , f m ) = K(h) if and only if h is a GCRC o/f 1; . . . ,f m . 

Proof. This proposition is proved in the univariate case in pQ but the proof can be 
extended to the multivariate case in a straightforward way. □ 

Proof of Proposition HH Firstly, we suppose that there exists a Liiroth's generator 
h = h\/h,2, where hi/h 2 is reduced. Then, by Proposition l28l / = uoh and g — voh 
where u, v G K(T) do not have a common right component of degree greater than 
one. Thus K(u(T),v(T)) = K(T). Then there exist Q U Q 2 6 K[U, V] such that 
Q 1 (u(T),v(T))/Q 2 (u(T),v(T))=T. 
Furthermore by Lemma [8l 

F a = f2{a)h(X) - h{a)h{X) = efj (h^X) - UteQQ) 

i 

where e G K and ti are the roots of 

/ 2 (a)ui(T) - fi(a)u 2 (T) =: u a , 

and 

G a = 92(a) 9l (X) - gi(a)g 2 (X) = e'J] (h^X) - Si h 2 (X)) 

i 

where e' G K and are the roots of 

g 2 {a)vi(T) - gi{g,)v 2 (T) =: v a . 

We get: h(a) is a common root of u a and v a . Thus hi (X) — h(a)h 2 {X_) divides F a 
and Gq. 

If f2(o L ).g 2 (g L ).Q2(u(h(a L )),v(h(a L ))^ ^ then /i(a) is the unique common root of u a 
and v a . Indeed if there exists another root x such that u a (x) — v a (x) — 0, then 
u(h(a)) = /i(a)// 2 (a) = u(x) and v(h(a)) = gi(a)/g 2 (a) = v(x). 
It follows: 

Qi(u(h(g,)),v{h(a))) Q 1 (u(x),v(x)) 

tl\Gi} — — = — = X. 

Q 2 {u{h{a)),v{h{g L ))) Q 2 (u(x),v(x)) 

Now we remark that if t ^ s then gcd(/ii + th 2 , hi + sh 2 ) is constant. 
We get then: gcd(F a , G a ) = hi(X) - h(a)h 2 (X). 
In the same way: gcd(Fb, G&) = hi(X) — h{h)h 2 {X). 

If h(a) ^ h(b), this gives the desired result, because K(h) = K(H) when H = U o h 
with U = (T — h{a))/(T - h(b)). 

Secondly, we suppose that there does not exist a Liiroth's generator. 
Then we have / = u o h and g — v o H , with h,H G non-composite and 

algebraically independent. 

Thus F a (X) = e.Ili (hi(X)-tih 2 (X)) as before, with h x {X) -Uh 2 {X) absolutely 
irreducible if ti a (hi, h 2 ). The condition ti $ cr(hi, h 2 ) means 

R(a) = Re ST (f 2 (a)ui(T) - fi(a)u 2 (T), S V ect hiM (T)) ± 0, 
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where Spect fcliha (T) = Uxe<r(h 1 ,h a )rK.( T ~ A )- 

In the same way, we have G a = e'. jX {Hi (X) — SiH 2 (2£)) with 

-Hi (20 — Si H 2 (X_) absolutely irreducible if 

S(a) = Res T {92{a)v 1 (T) - gi(a)v 2 (T), Speet HuIIa (T)) ? 0. 

Thus F a and G a have a non trivial common divisor if and only if there exist U, Sj 
and a E K \ {0} such that: 

(*) a(M20 - tih 2 (X)) = H^X) - s H 2 {X). 

In the same way, Fb and Gb have a non trivial common divisor if and only if there 
exists t't, s'j and a' S K\ {0} such that: 

(**) a'(h!(X) - t'-MX)) = HriX) - s'^X). 

(★) and (irk) give: 

a -atA (hx\ = (l -aA (E x 
a' -a't'J \h 2 ) \1 -s'J {H 2 , 

If 

D(a,b) = Res T (g 2 {a) Vl (T) - gi (a)v 2 (T), g 2 (b) Vl {T) - 9l (b)v 2 (T)) ^ 

then Sj ^ s'j and the previous system gives H = u o h, with degu = 1. Thus h and 
H are algebraically dependent and this is absurd. Thus F a and G a (resp. Fb and 
Gb) have no common divisor. 

Hence, if no Liiroth's generator exists and f 2 (a).g 2 (b)M(a).S(a)M(b).S(b).D(a,b) 
is not equal to zero, then gcd(F a , G a ) is constant and gcd(i^b, Gb) is constant. Thus 
the algorithm returns "No Liiroth's generator" . □ 

Remark 29. With the notations of the previous proof, we remark that a and b must 
avoid the roots of: f 2 (X), g 2 (X), h 2 (X), Q*{} : Q0) , 5(20) , ^(20, SQQ, and (a, b) 
must avoid the roots of hi(A)h 2 (B) - hi(B_)h 2 (A) and D(A,B). 
We can easily bound the degree of each polynomial: deg fa < d, deg gt < d, deg hi < 
d/2, degQ 2 < d(d-l) see Proposition 2.1], degi? < d(d 2 -l), degS" < d(d 2 -l), 
and deg D < d 3 . 

Then if K is "big enough" the open Zariski set U is not the empty set. 

Remark 30. In the algorithm Sederberg Generalized we cannot consider two random 
linear combinations of /i, f 2 and g\, g 2 . Indeed, with random linear combinations 
and with the notations of the previous proof, u a and v a do not have a unique 
common root in K. Thus with random linear combinations the strategy used in 
Proposition [25] is not valid. 

Proposition 31. 7/K is afield with at least (Ad + 2)d elements then the algorithm 
Sederberg Generalized uses 0(d n ) arithmetic operations. 

Proof. The computations of fi(a), <?i(a), fi(b), gi(l>) needs 0(d n ) arithmetic oper- 
ations. The complexity of an n-variate gcd computation needs (D(d n ) arithmetic 
operations. Indeed, as K is a field with at least (4c? + 2)d elements with Lemma 
6.44 in we can generalize to n variables the algorithm 6.36 presented in [3"9"] 
and obtain a result like Corollary 11.9 in [39 . This gives the desired result. □ 
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Remark 32. When it is possible, a polynomial generator is desirable. The algorithm 
Sederberg Generalized always returns a rational generator. We can test if we have a 
polynomial generator in the following way: We test if there exist a, (3 £ K such that 
aH a + f3 = Ht>. If such constants exist then H a (or Hb) is a polynomial generator. 
This improvement is correct because we have seen during the proof of Proposition 
[26] that H a = h% — h{a)h% and Hb — hi — h(b)h,2- Thus if a polynomial generator 
hi exists we have H a = hi — hi(a) and Hb = hi — hi(b). As gcd are known up to a 
multiplicative constant there exist a, /3 G K such that aH a + f3 = Hb- Conversely, 
if we have aH a + f3 = Hb then H a /Hb = u o H a with u — T/(aT + /?), thus 
K(H a /H b ) = K(H a ). 

The computation of a and /3 needs 0(d n ) arithmetic operations. Indeed, we solve 
a linear system with 0{d n ) equations and two unknowns. Thus we can find a poly- 
nomial generator with the algorithm Sederberg Generalized with 0(d n ) arithmetic 
operations. 

4.2. Another strategy based on decomposition. Now, we give another algo- 
rithm to compute a Liiroth's generator. Here we use the relation between decom- 
position and computation of a Liiroth's generator. 

Liiroth with Decomp 

Input: f(X) = fi/f 2 (X), g(X) = gi/g 2 (X) e K(X U . . . , X n ) two reduced rational 
functions, z := {a,b) <E K. 2n . 

Output: h(X) G KQC) such that K(f,g) = K(h), if h exists. 

(1) Decompose / with the algorithm Decomp, then / = u o h. 

(2) Compute v (if it exists) such that g — v o h. 

(3) If v do not exist then Return "No Liiroth's generator" , else go to [4] 

(4) Compute w the GCRC of u and v with Sederberg's algorithm. 

(5) Return w o h. 

Proposition 33. The algorithm Liiroth's with Decomp is correct for z satisfying 
the hypothesis of Theorem^ 

Proof. This algorithm computes a GCRC of / and g, thus by Proposition [251 this 
gives the desired result. □ 

Proposition 34. Under hypotheses (C) and (H), the algorithm Liiroth's with De- 
comp performs one factorization of a univariate polynomial of degree d over K plus 
a number of operations in K belonging to 0{d n ) if n > 3 or to 0(d 3 ) if n = 2. 

Proof. The first step of the algorithm performs one factorization of a univariate 
polynomial of degree d over K plus a number of operations in K belonging to <D(d n ) 
if n > 3 or to 6(d 3 ) if n = 2 by Proposition ITl 

With the strategy presented in Section [2Tl the second step can be done with <D(d n ) 
arithmetic operations. 

The last step can be done in an efficient probabilistic way, see [32]. The algorithm 
presented in |32) computes only two gcd's of univariate polynomials of degree lower 
than d. 

Then the total cost of the algorithm belongs to 0(d n ) if n > 3 or to 0(d 3 ) if 
ra= 2. □ 

Remark 35. During the algorithm Liiroth with Decomp we have to avoid the roots 
of nonzero polynomials considered in Remark [T3] and Remark [T5] because we use 
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the algorithm Decomp. Furthermore during the algorithm Luroth with Decomp, we 
use Sederberg's algorithm, this algorithm is also probabilistic and has in input two 
parameters x\,X2 G K. If x\ and X2 are not the roots of a nonzero polynomials 
then the output is correct, see [32] . 

Thus the nonzero polynomials are just the ones used for the algorithm Decomp and 
for Sederberg's algorithm. 

4.3. Computation of a Luroth's generator. 
Liiroth's generator 

Input: fi (X), . . . , fm.(X) e K(X), m reduced rational functions, 
z := z % , . . . ,z m G K 2n , n > 2. 

Output: h{X) £ K(X) such that K(fi, . . . , f m ) = K(h), if h exists. 

(1) Compute a Luroth's generator of K(fi,f2) with Sederberg Generalized ap- 
plied to fi, f2, with zi . 

(2) If a Luroth's generator h is found then go to step[3]else Return "No Luroth's 
generator" . 

(3) For i = 3, . . . , m, 

(a) Compute a Luroth's generator of K(h,U) with Sederberg Generalized 
applied to h, fj, with z\. 

(b) If a Luroth's generator H is found then h := H else Return "No 
Luroth's generator" . 

(4) Return h. 

Proposition 36. The algorithm Luroth's generator is correct for z satisfying the 
hypothesis of Theorem^ 

Proof. We just have to remark that K(fi, . . . , f<_i, fj) = K(fi, . . . , fj_i)(fj). □ 

Proposition 37. IfK has at least (Ad + 2)d elements, then the algorithm Luroth's 
generator can be performed with 0{md n ) arithmetic operations in K. 

Proof. We use m times the algorithm Sederberg Generalized. Thus, thanks to Propo- 
sition [3T] we get the desired complexity. □ 

Remark 38. During the algorithm Luroth's generator we can use the algorithm 
Luroth with Decomp instead of Sederberg Generalized. In the bivariate case, the 
complexity becomes then 0(d 3 ). In this case the algorithm is not softly optimal, 
but the algorithm can also return u such that / = u o h. 

We conclude that Proposition l36l and Proposition l37l prove Theorem 01 

5. Study of the Gutierez-Rubio-Sevilla's algorithm 

In this section we study the complexity of the decomposition algorithm given in 
[15] . More precisely, we explain how to modify it in order to get a polynomial time 
algorithm instead of an exponential time algorithm. 

5.1. Some preliminary results. The following lemma is a generalization of Lemma 

IE 

Lemma 39. Let h = hx/h^ be a rational function in KQQ, u — u\/u2 a rational 
function in K(T) and set f = u o h with f = fxjfi G KQQ. Let A, fj, G L, where L 
is a field and KcL. We have: 

M/i - A/2 = (pui - Xu 2 )(h).h2 CSU . 
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Proof. We have 

Vfi ~ A /2 _ u\{h) _ ^ u 2 (h) _ nui(h) - Xu 2 (h) 
h 112(h) u 2 (h) u 2 (h) 

Thus: (*) (ufa - \f 2 ).u 2 {h) = (jiui - Xu 2 )(h).f 2 . 
Furthermore 

h Mh) (Etohh\ht- i ).h^ 
where Ul (T) = ^ \ M T ) = Eto W*- 

Then f 2 = ( J2ilo bih\h 2 2 ~' l ).h™ ax ^ dl ~ d2 ' ^ because / is reduced and the degree of 
the right term of (**) is lower or equal to deg(/). 

It follows f 2 = u 2 (h).h™ & ^ dl d2 -°^ +d ' 2 — u 2 (h).h 2 csu , then thanks to (*) we deduce 
the desired result. □ 

Proposition 40. Let f E K.(X) be a rational function such that f = u o h and 
f = u o ip ! where u is a rational function in K(T), h a non- composite rational 
function and tp a rational function. 

Then tp is non-composite and there exists w G K(T) such that h — w o tp and 
degw = 1. 

Remark 41. w is not necessarily the identity. For example if u — x 2 + 1/x 2 and 
w = 1/x then uow — u. Thus we can get / = {uow)otp = uotp and / = uo(wotp) = 
u o h. See |16j for more statements on the particular situation u o w — u. 

Proof. We set u = u\ju 2 and tp = tpi/ip 2 . 

Let A,/i £l such that deg(/iiti — AM2) = deg it, by Lemma [39] we have 

dcg u 

M/l - A /2 = e JJ (hi - x { h 2 ), 
»=i 

where e G K and Xi G K are the roots of /xui — Xu 2 . 

We can suppose that hi — Xih 2 are absolutely irreducible and Xi ^ Xj if i ^ j. 
Indeed, the "bad" values of (fi : A) are (u 2 (x) : Ui(x)) where x G a(hi,h 2 ) and 
are the roots of R(fi, A) = Res([iui — Xu 2 , fj,u[ — Xu' 2 ). As cr(/ii, h 2 ) is finite and K 
infinite, we deduce that "good" values of (fj, : A) exist. 

We can also suppose that degt^i — Xitp 2 = degtp, because we just have to avoid a 

finite number of x%. 

Then Lemma [39] also implies 

dcg u 

M/i - A/ 2 = e J| (cpi - Xitp 2 ). 
i=i 

We have tpi — Xitp 2 is absolutely irreducible, else /i/i — A/2 has more than degu 
absolute irreducible factors: this is a contradiction with hi — Xih 2 being absolutely 
irreducible. 

Then tp is non-composite by Proposition [7] 

Furthermore, there exist ik,jk, with k = 1, . . . ,degu such that hi — Xi k h 2 equal 
tpi — Xj k tp 2 up to a multiplicative constant. As in the proof of Proposition 1261 it 
follows tp = w o h with w G K(T) and deg it) = 1. As h and tp belongs to KQQ we 
have w G K(T). (Indeed we just have to solve a linear system in K to get w.) □ 
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5.2. Study of the absolute irreducible factors of near-separated polynomi- 
als. The decomposition algorithm given in [15] is based on the following theorem; 
see [30) . In this subsection we improve this result. 

Theorem 42. Let f = f l / f 2 G K(X). 

f = uoh, with h = hi/h 2 if and only if H(X,Y) = hi(X)h 2 (Y) - h 2 (X)h 1 (Y) 
divides F{X,Y) = h(X)f 2 (Y) - f 2 (JQfi(Y). 

In the following we use a result due to Schinzcl. 

Definition 43. A rational function is reducible over K if the numerator in its 
reduced form is reducible over K. 

Lemma 44. Let ^(T,Y_) and f(X) be non-constant rational functions over K, the 
former of non-negative degree with respect to T and to at least one Yi . 
If the function 

is reducible overK. then f = uoh, u G K(T), h G 1KQ0 and ip(u(T),Yj is reducible 
over K. 

Proof. See [Ml Lemma 1]. □ 

Proposition 45. Let f — /1//2 G K(JQ, / = /1//2 G K(H) fee two non-constant 
rational functions. 

If f and f are non-composite then F(X_,Y_) — fi (X)h.(Y) — f? (X) fi (Y) is irre- 
ducible in K[X, Y]. 

Proof. We set ip(T, Y) = f(Y) - T. 
Then 

V ; / 2 (X)/ 2 (F) 

If we suppose -F(A, Y) reducible then f = u o h and ^(uf^)!!!) is reducible by 
Lemma [44] 

As / is non-composite deg u = 1 thus we can set m(T) = (aT + 6)/ (aT + /3). Then 
ip(u(T),Y_) is reducible means /i(Y)(a;T + j3) — f 2 (Y)(aT + b) is reducible over 
K. By Proposition [7] this is absurd because / is non-composite. Hence F(X_,Y2) is 
irreducible. □ 

Now we can improve Theorem 1421 

Theorem 46. Let f = /1//2 G K.(X) a non-constant rational function. 
If f = uoh, where u = u\/u 2 G K(T) and h = h\/h 2 G KQQ are rational functions, 
with degu > 2 and h non-composite, then the irreducible factors with the smallest 
degree relatively to X_ of 

F(X,Y) = h(X)f 2 {Y) - f 2 (X)h(Y) 

are of the kind 

H(X,Y) = hi[X)<pi,2(Y) - h 2 {X)<pi,i(Y), 

where tpi — tfi,i/ipi, 2 are non-composite rational functions such that h = w o ip^ with 
degw; = 1. 

Theorem [5] is a direct consequence of Theorem 1461 
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Proof. By Lemma l39l we have 

(*) f(x,y) = u fuf2 (h(x)).h 2 (2Q desu , 

where 

U flj2 (T) = fzQOMT) - h{Y)u 2 {T). 
As / = u o h, h(Y) is a root of Uf 1 j 2 . Then 

UfufAT) = (h 2 (Y_)T — hi(Y_))A(Y_, T), 

where A(Y,T) G K\Y, T] . Thus (*) implies hx(X)h 2 (Y) - h 2 (X)hi(Y) divides 
F(X,Y). 

Now, we suppose that tp(Y_) £ is another root of Uf x j 2 (T). Then 

«(<K£) = /Q0 = «(/»£))■ 

Thus, by Lemma l40l we have tp is non-composite and h — wotp with deg w = 1. As 
before, we can write U fuh = (tp 2 (Y)T -(pi(Y)).B(Y,T), where B(Y,T) eK[Y,T]. 
Thus ip 2 (Y)hi(X) - <pi(Y)h 2 (X) divides F(X,Y) by (*). 
Now, we write 

(**) U flj2 (T) = (^, 2 (Y)T - <p itl QO)- II C?QL,T), 
iei je.J 

where tp^ — <p% i/ipi 2 QD 1S a reduced non-composite rational function as explained 
above and Cj(Y,T) G K[Y,T] is irreducible with deg T C 3 > 2. 
We evaluate T to h in (★*) and multiply the result by h Acgu : 

u fuh {h{X))MK) dc&u = n(^(mQ0-^,i(r>2Q0) 

x ( n & h CK) )) ./i 2 (z) Ejej ej dogT Cj ■ 

The factors ipj ^(y)fti (X) — Pi,i(YJh 2 (X) are irreducible by Proposition H51 Fur- 
thermore, by Lemma l44l as /i is non-composite and Cj{Y_,T) is irreducible, we have 
C (^h(X)).h 2 (X_) dc&TCl is irreducible in K[X,Y}. 
We also have 

deg x Cj(Y,h(X))h 2 (X) de ZT c > = deg T C r degh 

> 2degh 

> deg x <Pi, 2 {Y)h x (X) ~ Vi A (Y)h 2 (X). 

Then H(X,Y) — (Pi ?(Y)hi (X) — i (F)fe^(X) are the factors with the smallest 
degree relatively to X_. □ 

5.3. Improvement of the GRS algorithm. Now we describe the decomposition 
algorithm presented in [T5] . 

GRS decomposition algorithm 

Input: f(X) = fi/f 2 (X), n>2. 

Output: u G K(T), h(X) G K(X) such that / — u o h, or "/ is non-composite". 

(1) Factor Let D = {Hi, . . .,H m } be the set of factors of F (up to 

product by constants). We set i = 1. 
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(2) If Hi can be written Hi(X,Y) = hi(X)h 2 (Y) - h 1 {Y)h 2 {X) then h 1 /h 2 is 
a right component for /. Then compute u by solving a linear system and 
Return u, h. 

(3) If i < m then i := i + 1 and go to step [2 else Return "/ is non-composite" . 
This algorithm has an exponential time complexity. Indeed, the set D contains 

at most 2 d polynomials, where d is the degree of /. 

However, we can improve this algorithm. Thanks to Proposition I45[ we remark 
that if / is non-composite then F is irreducible. Furthermore, if / = u o h with h 
non-composite, then H(X, Y) = fei (X)h? (Y) — hi (X)h 2 (Y) is an irreducible factor 
of F(X_,Y), by Theorem l46l Thus we have to study at most degF irreducible 
factors. Thus we can substitute the set D by the set of irreducible factors. (We can 
also substitute the set D by the set of irreducible factors with the smallest degree 
relatively to X). As Step [T] and Step [5] can be done in a polynomial time, it follows: 

Proposition 47. If in the GRS decomposition algorithm we set: "D is the set 
of irreducible factors of F ", then this modified algorithm has a polynomial time 
complexity. 

Remark 48. The bottleneck of this modified algorithm is the factorization of F. 
If we apply the deterministic algorithm proposed in [35] then the modified GRS 
decomposition algorithm uses 0(d 2n+u> ~ r ) arithmetic operations, where d is the 
degree of / and n the number of variables. 

Exemple 49. Now, we illustrate the GRS decomposition algorithm with / = u o h, 
where u — (T 2 + 1)/T, h = h 1 /h 2l and h x = X\ + X\ + 1, h 2 = iX 1 X 2 . h is a 
non-composite rational function. 

In this situation, we have the following factorization of F(Xi, X 2 , Yi, Y 2 ): 
F(X 1 ,X 2 ,Y 1 ,Y 2 )=3.H 1 (X 1 ,X 2 ,Y 1 ,Y 2 ).H 2 (X 1 ,X 2 ,Y 1 .Y 2 ), where 

H^XuXi^Yi) = XlY 1 Y 2 +XlY l Y 2 +Y l Y 2 -Y?XiX 2 -Y$X 1 X 2 -X 1 X 2 

= /n(x 1 ,x 2 )/ l2 (y 1 ,y 2 )-/i 1 (y l! r 2 )/ l2 (x 1 ,x 2 ), 

H 2 (X 1 ,X 2 ,Y U Y 2 ) = 1 + Xl + Xl + Y? + Yi + XfY? + X?Y 2 3 + X*Y? + X^ 
— 9Xi X 2 Y\ Y 2 

= ^ 1 (x 1 ,x 2 )^ 1 (y 1 ,y 2 )-/i 2 (r 1 ,y 2 )fo 2 (x 1 ,x 2 ). 

Then we can recover the decomposition / = u o h with the GRS decomposition 
algorithm. 
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